A route to computational chaos revisited: 
noninvertibility and the breakup of an invariant circle 

Christos E. Frouzakis 
Combustion Research Laboratory, Paul Scherrer Institute 
CH-5232, Villigen, Switzerland 

loannis G. Kevrekidis * 
Department of Chemical Engineering 
and Program in Applied and Computational Mathematics, 
Princeton University, Princeton, NJ 08544 

Bruce B. Peckham 
Department of Mathematics and Statistics 
University of Minnesota Duluth, Duluth, MN 55812 

February 1, 2008 

Abstract 

In a one-parameter study of a noninvertible family of maps of the plane aris- 
ing in the context of a numerical integration scheme, Lorenz studied a sequence of 
transitions from an attracting fixed point to "computational chaos." As part of the 
transition sequence, he proposed the following as a possible scenario for the breakup 
of an invariant circle: the invariant circle develops regions of increasingly sharper cur- 
vature until at a critical parameter value it develops cusps; beyond this parameter 
value, the invariant circle fails to persist, and the system exhibits chaotic behavior 
on an invariant set with loops [1^. We investigate this problem in more detail and 
show that the invariant circle is actually destroyed in a global bifurcation before it 
has a chance to develop cusps. Instead, the global unstable manifolds of saddle-type 
periodic points are the objects which develop cusps and subsequently "loops" or "an- 
tennae." The one-parameter study is better understood when embedded in the full 
two-parameter space and viewed in the context of the two-parameter Arnold horn 
structure. Certain elements of the interplay of noninvertibility with this structure, 
the associated invariant circles, periodic points and global bifurcations are examined. 
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1 Introduction 

In an insightful 1989 paper entitled "Computational Chaos: a prelude to computational 
instability" E. N. Lorenz [16j reported on a one-parameter computational study of the 
dynamics of the noninvertible map 



which arises from a simple forward Euler integration scheme (r being the time step of the 
integration) of the two coupled nonlinear ODEs 



These ODEs are obtained by starting with the familiar Lorenz system [T^: x = —a{x — 
y),y = —xz + px — y, z = xy — Pz, letting a — > cxd, and rescaling the variables and the 
remaining parameters. Lorenz fixed the value of a at 0.36, and varied the time step r. Al- 
though the corresponding differential equations (eq. exhibit an attracting equilibrium 
point, computer simulations indicated that the discrete approximation (eq. (^) progressed 
from exhibiting an attracting fixed point, to an attracting invariant circle (IC), to a chaotic 
attractor (termed computational chaos) as the time step r was gradually increased. Se- 
quences of bifurcations similar to those described by Lorenz have been observed for other 
noninvertible maps of the plane suggesting that this might be a universal "non- 

invertible route to chaos." Some of the bifurcations along the route — local changes of 
stability, homoclinic and heteroclinic tangencies, crises caused by an attractor interacting 
with its basin boundary — are also observed in families of invertible maps. New bifur- 
cations, however, unique to noninvertible families, as well as "invertible bifurcations with 
noninvertible complications," are also observed along the route. 

In this paper we perform a more detailed numerical study of Lorenz's family, but focus 
primarily on a narrow range of parameters which includes the breakup of the invariant circle 
and the identification of tcd (the greatest lower bound on the values of the parameter 
r for which the corresponding maps have "chaotic dynamics"), and Tchaos, the greatest 
lower bound on the values of the parameter r for which the corresponding maps exhibit 
a "chaotic attractor." Lorenz was interested in Tchaos (labelled Tfe in |Tni) as an indication 
of how poorly the Euler map L(^o.36,t) approximated the original differential equation. He 
observed that, in the transition from smooth IC to chaotic attractor, the process began 
with the circle developing features with increasingly high curvature. There appeared to 
be a critical r value at which the "IC developed cusps." After this value, the "chaotic 
attractor" suggested by computer simulations appeared to have loops and a Cantor-like 
structure. It certainly was no longer a topological circle. He claimed (correctly) that if this 
scenario did in fact occur, then this critical r value had to be an upper bound on Tchaos- 

Our numerical investigations suggest that a smooth IC does not persist all the way 
to a "cusp" parameter. Instead, the IC is destroyed in a heteroclinic tangency initiating 




(1) 



dx/dt = ax — xy 
dy/dt = —y + x^. 



(2) 



2 



the crossing of a branch of the stable manifold and a branch of the unstable manifold of 
a periodic (here period-37) saddle point. (As for invertible maps, chaotic dynamics are 
apparently present during this crossing, so this tangency is an upper bound on tcd-) It 
is the unstable manifold that subsequently develops cusps at a critical parameter value we 
label Tcusp, and loops beyond that value. A second manifold crossing results in the apparent 
appearance of the chaotic attractor, with loops inherited from the unstable manifold. We 
conclude that r^usp is not a single isolated bifurcation value separating smooth IC attractors 
from chaotic attractors. 

Rather, Tcusp is strictly above tcd and strictly below Tchaos- AH three parameter values 
are part of the transition from a smooth IC attractor to a chaotic attractor. Because a 
cusp on an IC with an irrational rotation number would necessarily force a dense set of 
cusp points on the IC, it would make the existence of the IC itself unlikely. We therefore 
expect the existence of an entire transition interval of parameters, as opposed to a single 
bifurcation parameter, to be the generic scenario for the transition between smooth IC and 
chaotic attractor in the presence of noninvertibility. 

Mechanisms for the breakup of ICs are relatively well established for invertible maps 
of the plane (see for example |31 IHl 1211 )• the invertible scenario which appears to 
us to most closely match the noninvertible scenario of this paper, an IC is born in a 
Hopf (also called Neimark-Sacker) bifurcation, grows in size, coexists with a periodic orbit 
(after a saddle-node birth of the periodic orbit off the IC), is destroyed in a first global 
manifold crossing (also referred to in the literature as a crisis |2S] where the IC collides 
with its basin boundary), and is reconstituted after a second global manifold crossing. The 
periodic orbit which persists through the destruction and reconstitution of the IC typically 
switches from outside (inside) the IC before it is destroyed to inside (outside) the IC after it 
is reconstituted. Chaotic invariant sets necessarily exist only during the manifold crossings. 
In contrast, the breakup of the IC in the noninvertible scenario is part of a transition to 
a chaotic attractor. In particular, the "chaotic attractor with loops," which appears in 
this noninvertible case, is not a feature which appears in invertible maps. In fact, cusps 
and loops on iterates of smoothly embedded curves are not possible for smooth invertible 
maps, but are common features on (invariant) curves of noninvertible maps r[T]. Further, 
bifurcations such as the manifold crossings (crises) can have additional complications due 
to the presence of noninvertibility. For example, the unstable "manifolds" involved in the 
global crossings may have self intersections and cusps. Or stable and unstable "manifolds" 
may cross transversely at one homoclinic point, but fail to preserve transversality at other 
points along the homoclinic orbit. 

In any case, the transition mechanism for the maps in eq. is different from the 
invertible transitions, and the noninvertible nature of the map plays an important role: 
the IC — for parameter values for which it exists — interacts with (actually crosses) the 
locus on which the Jacobian of the linearized map becomes singular. This locus (and 
its images and preimages) is crucial in organizing the dynamics of noninvertible maps; it 
is termed "critical curve" in the pertinent literature [HI UHl EH H] and constitutes the 
generalization (in two dimensions) of the critical point in unimodal maps of the interval 
jHllZI. In particular, as explained in [TT], the angle of intersection between a saddle unstable 
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manifold and the critical curve underpins the transition of the local image of the manifold 
from being smooth and injective to nonsmooth and injective (with a cusp) to smooth but 
not injective (with loops). 

As for invertible Hopf bifurcations (even in the noninvertible setting the dynamics 
near a Hopf bifurcation are necessarily locally invertible), an understanding of Lorenz's 
a = 0.36 one-parameter family is possible only by viewing it in the two-parameter Arnold 
horn context jHl IH 123 12H1 • We do this with L(^a,T), using the second parameter (a) already 
provided in eq. (P). This leads to an examination of the internal Arnold resonance horn 
structure, to be compared to and contrasted with the invertible case. The interiors of the 
resonance horns have features which vary significantly from those in the horns of invertible 
maps (see preliminary results in ^3 ^|). Although our understanding of this internal 
structure is far from complete, examination of the partial picture still provides insight 
into the bifurcations observed in the Lorenz one-parameter cut. Developing computational 
tools to further investigate the internal structure is part of our ongoing research. 

In this paper the observations of the transition in jTH] are briefly summarized in Sec. El 
and then revisited in detail in Sec. El In Sec. El we place our one-parameter cut in con- 
text by examining the "resonance horns" in the larger two-parameter space. We discuss 
related noninvertible issues, including computational challenges, in Sec. El and state final 
conclusions in Sec. IHl 

2 Observations 

Figure n briefly summarizes the relevant observations in pHj. For a fixed at 0.36, the 
differential equation (eq. (j2I)) has an attracting equilibrium point at {x,y) = (0.6,0.36). 
Its basin of attraction is the right half plane. There is a symmetric attracting equilibrium 
for the left half plane at (—0.6, 0.36). We will deal only with the right half plane attractors 
in this paper. The corresponding maps /^(o.se.r) all have a fixed point at (x, y) = (0.6, 0.36). 
In a one-parameter cut with respect to r, the fixed point is stable for small r, but loses its 
stability via Hopf bifurcation at r = 1.38889. The resulting attractor is initially a smooth 
IC (fig. UK, T = 1.55) but then develops progressively sharper features (fig. [Dd, t = 1.775, 
where the circle is still smooth even though it appears to have cusps), then self-intersects, 
suggesting chaos (fig.^, r = 1.785), and eventually exhibits clear signs of chaotic behavior 
(fig. mi, T = 1.91). Lorenz supports his claim of chaos with the computation of a positive 
Lyapunov exponent for r = 1.91 [11]. The last two figures are unusual in that the attractor 
apparently shows self-intersections, a clear sign of noninvertibility. All four figures were 
created by following a single orbit, after dropping an initial transient, from an initial 
condition near the fixed point. We use the term attractor loosely to denote the object 
generated by the computer via such a simulation. More formal definitions appear in the 
next section. 

Lorenz proposed the following scenario as the simplest possible transition of the at- 
tractor from a smooth IC to a chaotic attractor. As a parameter is increased (r in this 
example) and the IC crosses successively farther over the critical curve (called "Jq" below). 
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Figure 1: Phase portraits of the map L for a = 0.36 and (a) r = 1.55, smooth IC (b) 
r = 1.775, regions of sharp curvature appear on the IC (c) r = 1.785, self-intersections 
appear (d) r = 1.91, "fully chaotic" attractor 

causing it to develop successively sharper features. Simultaneously, its rotation number 
changes (increasing with r in this case). By restricting to r values for which the IC is 
quasiperiodic (and avoiding periodic lockings, defined below in Sec. 13. 1|) . one obtains a 
sequence of r values, each corresponding to a map with a smooth quasiperiodic IC, and 
limiting to Tchaos, the lower limit on chaotic behavior. At Tchaos, the attractor, which may 
or may not still be a topological circle, would develop cusps. Beyond Tchaos, the attractor 
would be chaotic, or at least have parameter values accumulating to Tchaos from above for 
which the corresponding maps would exhibit chaos. We will now examine the implications 
of such a scenario. 

3 Observations Revisited 

By performing a more detailed numerical study, we show below that, while Lorenz's sce- 
nario is the correct basic description of the transition, the IC actually breaks apart before 
the "cusp parameter" value is reached. The transition from smooth IC to chaotic attractor 
thus occurs over a range of parameter values, rather than at a single critical parameter 
value. We emphasize that Lorenz did not claim that his scenario did happen, only that 
there could be no simpler scenario. Our study shows that the proposed "simple" scenario 
does not happen in this particular family, nor should it be expected in the general transition 
from smooth invariant circle to chaotic attractor in the presence of noninvertibility. 
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Before beginning the description of our numerical investigation, we recall some termi- 
nology and preliminary results. 

3.1 Noninvertible preliminaries 

Let / : 3?^ — > 3ft^ be a smooth {C°°) map. Let x be a period-g point for /. As for invertible 
maps, we define the stable manifold of x to be W^{x) = {y G 3?^ : f'"^{y) — > a; as — > oo}. 
Since the inverse map is not necessarily uniquely defined, we define the unstable manifold 
of X as W^{x) = G 3?^ : there exists a biinfinite orbit {i/j} with = f{yi),yo = 
y and y-kq — x as A; — ^ oo}. The use of the term "manifold" is an abuse of terminology 
for noninvertible maps since, while the loeal stable and unstable manifolds are guaranteed 
to be true manifolds, the global manifolds are not [201 IHOl CI]- (Some authors address 
this problem by using "set" instead of "manifold" [H HH [13 IHl E21-) By iteration of 
the local smooth unstable manifold, the global unstable manifold can be computed as for 
invertible maps. If the manifold (or any curve) ever crosses the critical set (Jo, defined 
below) tangent to the zero eigenvector of the map at the crossing point, its image can have 
a cusp [m [33] ■ More degenerate cases are also possible, but not considered here. The 
global unstable manifold can also have self intersections. Thus, it is globally a smoothly 
immersed submanifold, with smoothness violated only at forward images of these special 
crossing points. This is discussed further below. The global stable manifold, due to multiple 
preimages, can be disconnected, fail to be smooth, and even increase in dimension |3!^ ITT] . 
The multiple preimages lead to interesting basins of attraction (see Section [HI and j21[2S])- 

A set 5 C 3?^ is invariant if f{S) = S. When S is invariant and a topological circle we 
call it an invariant cirele (IC). An invariant circle on which there exists a periodic orbit is 
called a frequency locked circle or circle in resonance. (The terminology is borrowed from 
return maps of periodically forced systems.) Typically we see a single attracting and a 
single repelling periodic orbit on a frequency locked circle, but multiple orbits are certainly 
possible. A periodic orbit is also called a periodic locking or a locked solution, although we 
usually reserve the use of periodic locking to refer to a periodic orbit on an invariant circle. 

We have already indicated above that we are using the term attractor informally to 
indicate whatever is suggested via computer simulations. More formally, we define a com- 
pact set K to be an attractor block if f{K) C K° (the interior of K). Then we define A, 
the attracting set associated with attractor block K hj A = fX^=i /"(-^) ^ [SO]- In the 
cases we consider, including for r in the range depicted in fig. [T] we take the attractor 
block K to be an annulus containing the "attractor" in each picture: inner radius a small 
circle around the repelling fixed point O, and outer radius a (topological) circle which is 
big enough to contain the attractor, but still in the right half plane. The attracting set 
will always be taken to be associated with this annulus. 

Any invariant circles or periodic points in the attractor block will also be contained in 
the attracting set A. So will the unstable manifolds of the periodic points. The attractors 
displayed by computer simulation are generally (an approximation of) some subset of A. 
If A exhibits at least some approximate recurrence property, such as a dense orbit (for 
example, a quasiperiodic IC), or a numerical orbit which appears dense (for example, a 
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Figure 2: (a) Invariant circle, F, one of its preimages, F^^, and the critical curves Jq and 
Ji (r = 1.775). The image of a "curve" crossing Jo (at A and B) is tangent to Ji (at 
L{A) and L{B), respectively), (b) "weak chaotic ring" possessing loops (r = 1.785); (c) 
enlargement showing the loop on the attractor at L{C). 

locked circle with a high period locking), then the attractor might be a good approximation 
of all of A. In other cases, such as a locked circle with a low period locking, simulation 
would only display the sinks. Computation of the unstable manifolds of the saddles (crucial 
to our numerical investigation) often gives a better picture of A than does a simulation. 
Note that a locked circle is still an invariant circle, even though the whole circle is not the 
attractor. The circle is smooth by the unstable manifold theorem through the saddles, 
but as with invertible maps, its smoothness where two branches of the unstable manifold 
meet at the sinks depends on the eigenvalues Unlike invertible maps, as indicated 
above, the smoothness of an IC can also be violated by the development of cusps. 

We now define the critical curve Jq (termed "curve C" in ^H] and often LC-i for Ligne 
Critique in the recent noninvertible literature fn\ I24j) as the locus in phase space where 
the linearization of the map becomes singular. We call its image Ji ("curve D" and LC, 
respectively). The two sets Jq and Ji, as well as additional images and preimages, are 
known to be a key to understanding the dynamics of a noninvertible map. Consider, for 
example. Fig. |2k showing an enlargement of the IC, F for r = 1.775, with three important 
additions: Jq, Ji, and F^^, one additional first-rank preimage of the invariant circle F. 
The map /^(o.se, 1.775) (abbreviated as L) has either one or three inverses depending on the 
phase-plane point in consideration (the term "Zi — map" has been proposed for such 
maps, |2HI)- The invariant circle F has three first-rank preimages: itself, the curve F^^ 
shown in the figure, and a third preimage F3 ^, further away in phase space (not shown). 
The geometry of the map can be visualized by first folding the left side of the phase plane 
along Jq and onto the right side. F and F^^ should now coincide. Next rotate roughly 
90° so that Jo maps to Ji. The image of F lands exactly on itself (before the folding), 
sending A to L{A) and B to L{B) in the process. Note that L maps points that are in 
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the intersection of the two regions bounded by T and to the region bounded by the 
pieces of F and Ji which connect at L{A) and L{B). That is, some points inside F are 
mapped outside F. Similarly, points inside F^^, but not F are mapped from outside F to 
inside F. Contrast this with the property that interiors of invariant circles for (orientation 
preserving) invertible maps are invariant. 

Note also that F is tangent to Ji at L{A) and L{B). The "folding" along Jq, typ- 
ically creates a quadratic tangency along Ji at images of curves crossing Jq, unless the 
curve crossing Jq does so at an angle tangent to the zero eigenvector at that point. This 
phenomenon has been extensively discussed in ^I]. Our explanation here is based on the 
understanding reached in that paper. As a parameter (in this case r) is varied, the map 
L(o.36,t} (again abbreviated as L), the entire attracting set A, the curves Jq and Ji as well 
as the intersection points of A with Jq (e.g. A) and with Ji (e.g. L{A)) vary. In our case, 
between r values for figs. and b, the invariant curve develops a "fjord' which pushes 
across Jq, creating two new intersection points of A with Jq, and two new tangencies with 
Ji at the two image points. We call the lower intersection point of the fjord C (see fig.Eb). 
The points C and L{C) vary with r as well. When the tangent of the "curve" in A at the 
intersection point C becomes coincident with the null vector of the Jacobian of the map 
at C, the image of A (A itself) acquires a cusp touching Ji at L{C) ^T]. As the parameter 
value T is further varied, the cusp becomes again a quadratic tangency, but A acquires a 
loop (figs. Eb and c). Iteration of this "loop" gives rise to an infinity of such loops (called 
"antennae" in ^B]) on A, suggesting chaotic behavior. 

Further, Lorenz gives a heuristic justification of a simple criterion for an attractor (not 
formally defined, but assumed by us to be a compact invariant set with some sort of recur- 
rence) to be a chaotic attractor (an attractor exhibiting sensitivity to initial conditions): 
the existence of two distinct points on the attracator which map to the same image point 
(also on the attractor) ^H]- In this context, the tracing out of an object with a single orbit 
suggests recurrence, and existence of a loop implies the existence of two points which map 
to the same image point. Thus, if this criterion is correct, attractors with loops would 
necessarily be chaotic attractors. 

3.2 Numerical investigation 

We set out to examine the transition from smooth IC to chaotic attractor in finer detail 
in figures El to ^2 Before we embark on a detailed description, we should state that the 
real picture is yet more complex, involving similar phenomena on an even finer scale. For 
example, the same phenomena which cause the destruction of the "large" IC could also 
destroy the "small" IC's which we refer to below. Nevertheless, this sequence of figures 
includes the essential features of the transition: coexistence of a periodic orbit with a 
nearby IC, destruction of the IC through a global bifurcation, appearance of loops on an 
unstable manifold, and the reappearance of an attractor, this time chaotic with loops. It 
might be useful to look ahead to the two-parameter bifurcation diagram in fig. El to put 
the one-parameter cut described in the current section in context. 

We start from the smooth IC at r = 1.775 (fig. At r ^ 1.776243 a period-37 
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Figure 3: Coexisting attractors for r = 1.7765: the attracting period-37 foci (•), and the 
IC. The period-37 saddles (x) are also shown, along with both branches of the unstable 
manifold 5*5. One branch approaches the focus F^, while the other branch approaches the 
IC. 

saddle and node pair appears in a saddle-node (turning point) bifurcation. The bifurcation 
occurs off' of the IC, and results in the coexistence of two attractors: the stable invariant 
circle, F, and an attracting period-37 sink. Although the sinks have real eigenvalues at 
the saddle-node bifurcation, the eigenvalues very quickly become complex and they have 
been labelled in fig. Ob for r = 1.7765 by for "focus." The separatrix between these 
two coexisting attractors consists of the stable manifolds of the period-37 saddles. Global 
unstable manifolds have been systematically constructed as (37th) forward iterates of local 
unstable manifolds as for invertible maps. In this case, we have plotted only the two 
branches emanating from saddle S^. The right branch of the global unstable manifold 
directly approaches F^. This is clear in the blowup of fig. Eh- The left branch very quickly 
approaches the IC F. Successive 37th iterates travel slowly in a clockwise direction around 
the IC, eventually showing a good approximation to the whole IC. The approach to the 
IC is somewhat clearer in fig. Eb, where the branch appears to "join" the IC just to the 
left of S'5. The unstable manifold then continues almost "on" the IC: past F2S, S2S, and 
Fi4, evolving around the IC, reentering from the right, just below F5, before completing 
the (approximation to the) IC just to the left of S'5. 

Notice also that 18 of the period-37 saddle-node pairs are "inside" (those marked by 
numbers inside F) and the rest "outside" the invariant circle; this is another feature of 
noninvertible maps [llj). We discuss this again in Sec. 

As T increases further, the consequences of global bifurcations are observed; these bi- 
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Figure 4: During the "first" manifold crossing (r = 1.77688). a) tlie period-37 saddle 
and sink solutions have been plotted with one branch of the global unstable manifold 
(the branch involved in the crossing) which emanates from each of the 37 saddles, b) 
Both branches of the unstable manifold emanating from the saddle labelled Sig have been 
plotted (but not the branches emanating from 5*5, or any other saddle). 

furcations involve interactions between the stable and unstable manifolds of the saddle 
period-37 points. The first such global bifurcation is evident in the sequence of figures El 
through El the left branch of the unstable manifold of 5*19 misses to the left in fig. El and 
approaches the IC (the curve coming into the inset from the right and below is prac- 
tically coincident with this unstable manifold branch); it makes the "saddle connection" 
with S5 in fig. Hj it misses 5*5 to the right in fig. Eland approaches F5. 

The global crossing actually "starts" with a heteroclinic tangency (between a branch 
of a stable manifold of a period-37 saddle and and branch of an unstable manifold of 
another saddle on the same orbit) which destroys the IC at r 1.776878. This tangency 
is also called a crisis since as r increases toward the tangency parameter value, the 
IC (approached by the unstable manifold) and its basin boundary (the stable manifold) 
approach each other, colliding at the tangency parameter value. The attracting set during 
the crossing appears — as for invertible maps — to be a complicated topological object 
which at least includes the closure of the unstable manifolds of the period-37 saddles. 
Numerically, most orbits appear to approach the period-37 sink (fig. EJ. Thirty-seventh 
iterates of each of the 37 branches of the unstable manifold plotted in fig. evolve clockwise 
toward the "next" saddle, suggesting that we are at least close to a heteroclinic crossing 
of stable and unstable manifolds. Although the stable manifolds of the saddles cannot be 
realistically computed because of the exponential explosion ( "arborescence" ) of the number 
of period-37 preimages, the crossing of the stable manifold is evident in the inset (fig. Eb). 
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The left branch of the unstable manifold of 5*19 heads toward saddle 5*5, and although 
it is not clear from the figure, it begins to alternate back and forth, with folds on the 
right approaching F5, and folds on the left approaching 5*28- The portion approaching 
5*28 folds again, with one side heading toward F28 (visible) and the other side heading 
toward S14 (only the start is visible). Further iterates (not shown) would result in this 
one branch of the unstable manifold of S'19 eventually evolving all the way around the 
attracting set suggested in fig. H^, repeatedly visiting every saddle and every node. As 
with invertible maps, this allows for long chaotic transients, if one could appropriately 
choose initial conditions on the unstable manifold. Most orbits, however, eventually are 
attracted to the period-37 sink. 

We note that a transverse crossing of stable and unstable manifolds in the noninvertible 
setting does not imply transverse crossings along the rest of the orbit. For this reason, 
homoclinic or heteroclinic crossings do not necessarily imply the existence of a "horseshoe" 
or the chaotic transients referred to above. However, Sander (HH] shows that the failure 
of a transverse crossing to provide a horseshoe is a co dimension-two phenomenon, so we 
assume that our transverse crossings in our one-parameter study avoid such a problematic 
point, and behave as in the invertible case. Thus we expect the parameter value at the 
tangency initiating this heteroclinic crossing to be an upper bound on tcd, marking the 
first appearance of chaotic dynamics. The invariant set with the chaotic dynamics does not 
appear to be an attractor, however, so the chaos is numerically evident only as transient 
behavior. 

After the crossing, which ends in a heteroclinic tangency aX r ^ 1.776881, both branches 
of the unstable manifolds of the period-37 saddles approach the period-37 sink. The only 
attractor that remains is the stable period-37 sink, shown along with the period-37 saddles 
in fig. Elfor r = 1.1778. The attracting set, however, is a "new" invariant circle, with a 
period-37 locking, formed by connecting the saddles and nodes with the unstable manifolds 
of the saddles. This invariant circle is seen in the blowup of fig. Eb- It is smooth except 
at the sinks, where two different branches connnect and the smoothness depends on the 
eigenvalues jlj. Here the IC is only since the sinks are foci, with complex eigenvalues. 
There are apparently no cusps on the unstable manifolds at this parameter value. 

As T is increased, a secondary Hopf bifurcation of the period-37 foci at r = 1.779444 
renders them unstable. This can be seen in fig. IHlfor r = 1.78, where the unstable foci, 
are surrounded by stable period-37 invariant circles. It is important (and barely visible 
in the figure) that the unstable manifold intersects the small invariant circles (and itself!) 
repeatedly. This is a reminder of the noninvertibility of the map, and our acknowledged 
abuse of the term "manifold." See [llj for a discussion of this type of unstable manifold 
self- intersection. 

A number of high-period lockings on the period-37 invariant circles appear for very 
narrow ranges of the parameter. For example, at r = 1.7835, a period-259 saddle-node 
pair (formed after a period-7 locking of the period-37 invariant circles) is visible (fig. Ej); 
assuming the period-7 locking is on the small ICs, these small ICs are now frequency locked 
circles, composed of the chain of the unstable manifolds of the period-259 saddles (not 
shown). The period-37 unstable manifolds asymptotically approach the small ICs. The 
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Figure 5: After the first manifold crossing (r = 1.778): tlie invariant circle as an attractor 
has disappeared, but a new invariant frequency locked circle has formed as an attracting 
set. Both branches of the unstable manifold of the two adjacent period-37 saddles (5*5 and 
5*19) are plotted. The right branch of the unstable manifold of and the left branch of 
the unstable manifold of -S'lg meet at the attracting focus F5. 
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Figure 6: After a secondary Hopf bifurcation of the period-37 foci: period-37 invariant 
circles surround unstable foci (o) at r = 1.78. Both sides of the unstable manifold of the 
period-37 saddles — shown here only for — asymptote to the stable period-37 invariant 
circles. 
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Figure 7: Locking on the small IC (r = 1.7835). Both branches of the unstable manifold 
of 5*5 are plotted. Self-intersections are apparent on the right branch. Also shown are 
period-37 saddles (x), period-37 repelling foci (o), and, in (b), period-259 sinks (•), but 
not the period-259 saddles. 

right branch clearly intersects itself and the IC it approaches (fig. UJp)- Most individual 
orbits eventually approach the attracting period-259 sinks. 

A second heteroclinic manifold crossing involving the "other" side of the global unstable 
manifold of the period-37 saddle "starts" at r ^ 1.78373. The bifurcation is apparent in 
figures [7| through El the right branch of the unstable manifold of swings around F^, and 
misses Sig to the left in fig. [7J approaching the small IC around F^; it makes the "saddle 
connection" with 5*19 in fig. |H1 it misses S'lg to the right in fig. IHl with successive iterates 
following essentially "on" the attractor with loops. 

As for the earlier manifold crossing, understanding the saddle connection at r = 1.784 
requires some explanation. At this parameter value (fig. |H1), a period-8 resonance on the 
small invariant circles results in the period-296 points apparently being the only attrac- 
tors. The left branch of the unstable manifold of very quickly approaches the small 
IC surrounding F2S and continues counterclockwise around it; it is not a homoclinic loop, 
although it appears close to being one. The right branch heads for S'lg after looping clock- 
wise about 3/4 of the way around F5 (including the two high curvature "corners"); the 
global crossing is apparent since part of the manifold iterates to the left of S'lg around the 
IC (we iterated only half way around the IC), and part of the manifold iterates to the right 
of S'lg toward the right edge of the inset; this part to the right of S'lg continues the "long 
transient" behavior. It is important to notice that although both branches of the unstable 
manifold of S'lg are being followed by parts of the right branch of the unstable manifold of 
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Figure 8: During the "second" manifold crossing (r = 1.784). (a) Unstable manifolds of 
all 37 saddles are plotted, (b) Only the two branches of the unstable manifold of 5*5 are 
plotted; the right branch is the heteroclinic branch. For clarity, neither the period-296 
saddles nor foci locked on the small ICs are plotted. 

neither branch of the unstable manifold of 5*19 has been directly plotted in this figure. 

The slope of the unstable manifold at the point of intersection with Jq continues to 
change as r is increased. At r = Tcusp ~ 1.78428 the slope becomes equal to the slope of the 
eigenvector corresponding to the zero eigenvalue, and cusp points appear on the "manifold" 
(fig-Efc)- Beyond this "cusp parameter" value, the unstable manifold contains loops. Since 
'Tcusp occurs while the unstable manifold is still intersecting the stable manifold, the picture 
for r just beyond Tcusp, but not beyond the heteroclinic tangency which marks the end of 
the stable/unstable manifold crossing, is more complicated than might be expected in 
other families. In our case, the newly developed loops are damped out as they eventually 
approach the small period-37 ICs (possibly with locked solutions on the small ICs). Of 
course, the loops could follow a long transient before being attracted to the ICs. We do 
not have a figure corresponding to such a parameter value. 

With further increase of r, we pass the "other end" of the manifold crossing (beyond the 
final heteroclinic tangency) and an apparently chaotic attractor appears. Thus, this final 
heteroclinic tangency parameter value appears to be an upper bound on Tchaos, marking 
the first appearance of a chaotic attractor. At r = 1.785 the resulting attractor is seen 
in fig. inito coexist with a stable period-407 (born at a period-11 locking on the period- 
37 small invariant circles): the left branch of the unstable manifold of 6*5 approaches 
the nearby small IC with the period-407 locking; the right branch approaches the large 
amplitude "messy" attractor that self-intersects and exhibits the loops it has inherited 
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Figure 9: (a) Coexistence of period-407 attracting orbit with a "weakly chaotic ring" 
(r = 1.785). Plotted are the period-37 saddles (x), period-37 repelling foci (o), period-407 
attracting nodes (•), and the unstable manifolds of the period-37 saddles, (b) enlargement 
with only the two branches of the unstable manifold of 5*5 plotted, (c) Further enlargement 
of a portion of the saddle global unstable manifold in (b) at r = 1.785 (solid line) with 
loops; the analogous unstable manifold for the approximate "cusp parameter" : r = 1.78428 
(dashed line). Not shown are the period-245 saddles and nodes which exist on the small 
ICs at this cusp parameter value. 

from the global unstable manifold. Such an attractor has been called a "weakly chaotic 
ring" dm I24p. Note that Lorenz plotted this weakly chaotic ring attractor by iterating 
a single orbit starting near the fixed point (fig. ^); we plotted the same attractor by 
following a branch of the unstable manifold of the period-37 saddle orbit (fig. EI)- Tracking 
the periodic saddles and their manifolds also allowed us to see the coexisting small IC with 
the period-407 attractor which was not apparent in fig. 

The small amplitude invariant circles now begin to shrink, disappearing at r ^ 1.78626168 
in a second secondary Hopf bifurcation; the period-37 foci from return from sources to 
sinks; these continue to coexist with the large attractor that has become very "messy" 
with many more points of self-intersection (fig. llOb for r = 1.7866). The stable period-37 
foci approach the period-37 saddles in anticipation of the following period-37 saddle-node 
bifurcation, which happens at r ~ 1.7866262. 

Finally, at r = 1.788 the saddle-node bifurcation of the period-37 solutions has occurred, 
and the pair has disappeared. Note that (as with invertible maps) the pairings have 
switched. For example, disappears with F2S, while F5 disappears with Sig. Only the 
large chaotic attractor is left, seen in fig. ^2 

Recap. The transition from smooth IC attractor (r = 1.75, figs, ^p, to nonin- 
vertible "messy" chaotic attractor (r = 1.788, fig. [TT|) is not a result of a smooth change 
of a smooth, large amplitude invariant circle which develops cusps at a critical parameter 
value. Instead, the transition occurs via a sequence of bifurcations, notably involving the 
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Figure 10: Coexistence of period-37 attracting orbit with a "weakly chaotic ring" (r = 
1.7866). (a) Plotted are the period-37 saddles (x), period-37 attracting nodes (•), and the 
unstable manifolds of the period-37 saddles, (b) Enlargement with only the two branches 
of the unstable manifold of 5*5 plotted. The left branch immediately approaches the period- 
37 stable orbit (-F28); the right branch approaches the "big" attractor: it iterates off the 
figure to the right, follows counterclockwise around the attractor seen in (a), and returns 
in the bottom of (b), coming in just to the right of and S^; its loops are larger than in 
fig.Eb. 
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Figure 11: Only the chaotic attractor remains (r = 1.788). It was approximated by 
computing a single orbit. The period-37 fixed points have disappeared in a saddle-node 
bifurcation. 

interaction of the IC with nearby periodic orbits and the unstable manifolds of the saddles. 
Specifically, the IC as an attractor is destroyed in a global manifold crossing that leaves an 
attracting periodic orbit as the only attractor. The global unstable manifolds of the sad- 
dle period-37 points are the invariant objects that first develop cusps (at Tcusp ~ 1.78428, 
fig. int) and later loops (beyond Tcusp)- When the cusps and loops first appear, they "damp 
out" as they approach the attracting periodic orbit (or small IC). When the large ampli- 
tude attractor reappears at after the second global manifold crossing, it is approached by 
the branch of the unstable manifold which now has loops. These loops are inherited by 
the large amplitude attractor which is now apparently chaotic. 
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Figure 12: (a) The start of a two-parameter bifurcation diagram for the map L(^a,T), con- 
taining the fixed-point Hopf bifurcation curve HBi, and three pairs of saddle- node curves 
emanating from HBi corresponding to period-5, period-6, and period-37 resonance horns. 
The x's are labelled by the phase portraits in fig. ^ (b) An enlargement of the period-37 
resonance horn where it is traversed in our one-parameter cut. The two period-37 saddle- 
node curves (SN37) mark the upper and lower boundaries of the resonance horn. Two equal 
eigenvalue curves and a secondary Hopf curve (HB37) subdivide the horn. The parameter 
values marked with x's are labelled by our corresponding figure number (s). 

4 The larger picture. 

While a rich dynamical picture occurs along our one-parameter family with respect to r, 
it is clear that it only constitutes a part of a larger, actually two-parameter scenario. This 
is consistent with existing two-parameter descriptions of Hopf bifurcations for invertible 
maps [SI m El EEl EZl EH] • The primary features in the invertible case are resonance horns — 
inside which the corresponding maps have a periodic orbit of a certain period — emanating 
from the primary Hopf bifurcation curve. Since the fixed point at a Hopf bifurcation has 
eigenvalues on the unit circle, the IC close to its birth lies far away in phase space from Jq, 
where the map has a zero eigenvalue. Thus, we expect the two-parameter bifurcation dia- 
gram for the noninvertible Hopf bifurcation to strongly resemble the analogous bifurcation 
diagram for the invertible case — at least near the Hopf bifurcation curve. 

We have embedded our one-parameter family (a = 0.36 in eq. (^) in the two-parameter 
family i^(a,r)- Fig. IT^ shows part of the "big picture" in the r-a plane: the period-one 
Hopf bifurcation curve together with period-37, period-5, and period-6 resonance horns 
emanating from it. The bifurcation points we encountered along the one-parameter cut 
"continue" to bifurcation curves in the two-parameter setting. The curves displayed in this 
figure are only a tiny sampling of the dozens of bifurcation curves we have computed, which 
are, in turn, only a tiny sampling of the bifurcation curves that can be computed for this 
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family. (In particular, there are saddle-node curves in addition to those emanating from 
the primary Hopf bifurcation curve. See fI2\ 1^ for descriptions of such intraresonance 
horn features for invertible maps.) 

As we move higher up in the horns — above the "IC tangent to Jo" curve in fig. IT^ — 
the circle grows and starts to interact with Jq and thus becomes "truly" noninvertible. It 
is above this curve where the corresponding maps can exhibit noninvertible behavior, such 
as mapping some of the interior of an IC to its exterior and vice versa (recall fig. |2t). More 
"interesting" behavior, such as coexisting attractors or the breakup of the IC, is suggested 
when the resonance horns start to overlap. In our one-parameter cut this overlap happens 
at least by the first period-37 sadde-node bifurcation at r ~ 1.776243, although the three 
horns we have displayed do not overlap until higher r values. 

The enlargement in fig. IT^ shows three of the multitude of bifurcation curves which 
subdivide the period-37 horn into "inequivalent" regions: two "equal eigenvalue" curves 
(just inside the horn), marking the transition between real and complex eigenvalues for 
a period-37 orbit, and a secondary Hopf curve (HB37), marking the birth of the "period- 
37 small ICs." Somewhat unsatisfying in fig. IT^ is the fact that none of the displayed 
bifurcation curves is unique to noninvertible maps. We remark that we did compute an 
"eigenvalue zero" curve (for the existence of a period-37 point with an eigenvalue of zero), 
but we did not display it in fig. IT2b because it had no direct bearing on the bifurcation 
sequence we described (but see such a curve in a period-6 horn in fig. UHl). However, this 
curve certainly has an indirect bearing on the bifurcations, because it is related to the 
crossing of Jq by the attracting set (which includes any ICs, periodic orbits, and their 
unstable manifolds). If the attracting set never interacted with Jq, the bifurcations could 
be produced with invertible maps. 

Other bifurcation curves, which our current computational tools did not handle easily, 
would separate the other x's in fig. from each other. For example, the x labelled E] 
is in a thin heteroclinic region; heteroclinic tangency curves separate it from 01 below and 
El above. Similarly, the x's labelled |H1 and IHt are in a thin heteroclinic region, separating 
them from 13 below and ini above. A "cusp curve" (unique to noninvertible resonance horns) 
runs through IHt- The x labelled [3 is in a thin "period-259 secondary resonance horn," 
which emanates from the secondary Hopf curve HB37 and corresponds to a locking on the 
small period-37 ICs. 

A more complete knowledge of the internal structure of the "noninvertible" resonance 
horns in the two-parameter setting is clearly the key to extending our current understanding 
of the transition from a smooth invariant circle to a chaotic attractor. 

5 Discussion 

5.1 Typical one-parameter cuts 

A natural question which arises following the detailed numerical investigation of Sec. El is 
how "typical" are the features of that specific one-parameter family. We take the liberty 
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in the next few paragraphs of using the intuition we have gained via our numerical explo- 
rations, including the partial two-parameter bifurcation diagram of the previous section, 
to provide at least a partial answer to this question. Even though somewhat speculative, 
we feel that these observations are useful for placing the one-parameter cut in a broader 
context. 

A key observation accompanying our understanding of the IC breakup is that an IC 
with a cusp and an irrational rotation number would have the cusp iterate forward to a 
dense set of image points, which would also be cusps. Even if such an object is topologically 
possible, it appears to us to be highly unlikely. Because of this, we expect the invariant 
circle to break up before it reaches a parameter value where it develops a cusp. Thus, the 
IC must interact with some other object in the attracting set. In our experience, the other 
object in this "pre-cusp" breakup is a nearby periodic orbit and its accompanying stable 
and unstable manifolds. The related manifold crossings appear to be the "natural" way to 
first destroy an existing IC and to later reconstitute "itself" as a chaotic attractor. 

This observation is consistent with the sequence of transitions we described in Sec. El 
Extending this one-parameter transition interval to the two-parameter setting, we imagine 
a transition band, roughly parallel to the TC tangent to Jq' curve, and passing between pa- 
rameter values labelled 'lb' and 'Id' in fig|12b. Below the lower boundary of the transition 
band a smooth invariant circle exists (possibly locked, possibly coexisting with periodic 
orbit (s)). Above the upper boundary of the transition band, the corresponding maps are 
chaotic. Loosely, the transition band is a "thick version" of the "cusp curve" suggested by 
Lorenz. 

Further supporting this transition band scenario is the prominent "opening" of the 
resonance horns in fig. IT2k . This would be even more prominent if more resonance horns 
were plotted: most of the resonance horns appear to open up and overlap near the top of 
the upper boundary of the transition band. This is suggestive of, or at least consistent 
with, chaos. 

The existence of the transition band also suggests why the period-37 orbit might have 
been the specific orbit that was involved in our one-parameter scenario: that was the horn 
which was passed through as the transition to chaos band was crossed. Vertical parameter 
cuts for values of a other than a = 0.36 could pass through a different periodic horn 
during this narrow transitional parameter interval. We note that Lorenz specifically chose 
a = 0.36 because it appeared to "avoid periodic windows." In the two-parameter context, 
this can be interpreted as avoiding low-period resonance horns. Thus, we might expect a 
"typical" one-parameter cut to often interact with a periodic orbit of a period less than 
37. 

Alternative viable bifurcation sequences are suggested by the blowup of the period-37 
horn in fig. IT^ . For example, moving the cut to the left of a = 0.3599 could avoid the 
secondary Hopf bifurcation curve HB37, and thus avoid the small invariant circles that 
appeared in the a = 0.36 cut. Other differences in varying a would be more obvious if 
we were able to compute the region(s) of the period-37 horn corresponding to the global 
bifurcations, and curves corresponding to the development of cusps. (We expect to be able 
to do this in the future.) We know the x labelled '9c' in fig. IT^ is both on a cusp curve 
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Figure 13: Period-6 resonance horn and related bifurcation curves. 

(in the parameter space) and in a heterochnic region. Moving the vertical cut either left 
or right would probably prevent these two from coinciding. That is, the fact that the cusp 
parameter value occured during the manifold crossing in the a = 0.36 cut is probably not 
typical of most one-parameter cuts. If the cusp parameter value occured first (with respect 
to r), then the phase portrait for the first appearance of loops would be simpler: the loops 
would simply approach the small attracting ICs. (If we encountered the cusp parameter 
outside the secondary Hopf curve, it would be still simpler: the loops would damp out at 
the attracting period-37 nodes. See fig.HSlfor comparison.) On the other hand, if the cusp 
parameter value occured after the heterochnic region had been crossed, we claim (with no 
proof) that there must be another periodic orbit interacting with the IC as well. Perhaps 
its unstable manifold will be the one to develop a cusp and /ater undergo a global manifold 
crossing to reconstitute the chaotic attractor with loops. 

5.2 Related noninvertible features 

The relatively high period of the locking encountered along the one-parameter cut above 
makes the elucidation of the phase and parameter space structures cumbersome. To illus- 
trate some of these phenomena more succinctly we perform computations in the neighbor- 
hood of a simpler (period-6) resonance horn. 

Fig. ^1 is a blow-up of this period-6 resonance horn, with two additional types of 
bifurcation curves plotted inside: (a) two global bifurcation curves (approximating a thin 
heterochnic manifold crossing region associated with the stable and unstable manifolds 
of the period-6 saddle points) and (b) the curve on which the periodic points acquire an 
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Figure 14: Inside the horn: Period-6 saddle-node pairs inside and outside the invariant 
circle (a = 0.465, r = 1.472); the hatched region inside F is mapped to the grey region 
outside. 

eigenvalue zero (i.e. one periodic point lies on Jq and one on Ji) [121 Il2j. This latter 
curve provides an indication of the interaction of the dynamics with noninvertibility. An 
additional curve ('IC tangent to Jo'), extending beyond the resonance horn, is plotted. 
This is the locus, in parameter space, where the invariant circle starts interacting with 
noninvertibility, by first touching the critical curve Jq. Note that the 'eigenvalue zero' 
curve is tangent to the 'IC tangent to Jq curve. 

In fig. El we show a phase portrait for a parameter value inside the period-6 horn and 
above the 'IC tangent to Jq' curve ( (a, r) = (0.465, 1.472) ). The invariant circle crosses Jq 
transversely at two points A and 5; it is easy to see that it will then become quadratically 
tangent to Ji at the image points L{A) and L{B) (see fill I24p . Also illustrated in the 
figure is the relative location of the attracting invariant circle and one more of its three 
first-rank preimages. Clearly, one of the period-6 saddle-node pairs is located inside the 
attracting invariant circle while the remaining five are located outside. (Note that iterates 
of points to the left of Jq switch sides of the IC; thus, there must be an even number of 
periodic points to the left of Jq.) We have also illustrated the portion of the interior of 
the invariant circle that is mapped "outside" after one iteration of the map, as well as the 
portion of phase space outside the invariant circle that gets mapped inside. 

Our next observation is for another parameter value inside the period-6 horn: (a, r) = 
(0.4838, 1.409). In fig. ^1 we have an example of an unstable manifold of a period-6 saddle 
with a cusp which damps out at the period-6 attracting orbit. Recall that in our one- 
parameter cut, the cusp on the period-37 unstable manifold developed in the middle of 
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Figure 15: Simple damping out of a cusp at an attracting period-6 node ((a,r) = 
(0.4838, 1.409) ). (b) The upper right cusp maps (under L^g ^ggg ^ ^gg-,) to what looks like the 
end of a line segment emanating from A^^i. The whole unstable "manifold" has effectively 
collapsed onto this line. Successive (sixth) iterates of the cusp approach A^^i along the line. 

a heteroclinic crossing. Iterates of the cusp could therefore follow a long transient before 
eventually being damped out as it approached a period-37 sink. Fig. ^Jis a simpler version 
of the damping out of cusps on an unstable manifold in the approach to a sink. 

Finally, we have made an attempt to illustrate the complexity of the basins of attraction 
of the coexisting attractors: the stable period-6 points and the stable invariant circle. While 
for the purposes of this paper we will not dwell on the details of this complicated structure 
(for a discussion of the basic building blocks of such a picture see again [llj as well as 
[21 ini Cni UHl [I2l IHj), we emphasize a few key elements: 

(a) the basins of attraction are not simply connected; 

(b) the stable "manifolds" consist of a large number (possibly, even an infinity) of 
branches (resulting from the multiple backward in time trajectories). Only the com- 
puted "immediate" local manifold of the period-6 saddles is shown in fig. ^1 the 
images and preimages of them form the boundaries of the grey disconnected "is- 
lands" that belong to the basin of attraction of the period-6 foci. 

(c) notice that the stable set intersects Jq. The images and preimages of the intersection 
points have been called "points of alternance" ^H] on parts of the stable manifolds: 
these are points (on Jq and its preimages) where two preimages of a branch of the sta- 
ble set come together; the "direction of relative movement" of nearby points switches 
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Figure 16: Basin of attraction of the period-6 foci marked by filled circles (grey region) 
and of the coexisting invariant circle (the white complement) at (a, r) = (0.45, 1.563). The 
computed "immediate" stable manifolds of the period-6 saddles (marked by crosses) are 
also plotted. 

at such points of alternance. 
5.3 Computational issues 

As a final comment on the current state of our work on noninvertible systems, we note 
that computational tools for the accurate approximation and visualization of objects such 
as stable and unstable manifolds are crucial in both understanding what is possible and in 
analyzing specific examples. The possibility of multiple inverses can cause an exponential 
explosion in the number of distinct segments of an invariant curve or of stable manifolds 
of periodic points. The same is true for the critical set: all the images and preimages 
of the critical curve Jo- In three dimensions the computational issues are even further 
complicated: the Jq curve becomes a surface, and its iterates and preiterates must be con- 
structed. Standard one-dimensional continuation tools have to be extended (to simplicial 
continuation, or through software like PISCES |31]) to effectively compute, triangulate and 
visualize these surfaces and their intersections. Constructing scientific computing, database 
management and visualization tools for these objects in two and three dimensions is an im- 
portant task, and to our knowledge, an open area of research. Development of these tools 
would facilitate the analysis of many physical systems that are inherently noninvertible: 
discrete-time adaptation, control and decision problems, ranging from Model Reference 
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Adaptive Control systems to signal processing applications -in addition to more tradi- 
tional noninvertible models arising in population dynamics and economics. These tools 
will also be applicable in the case of multiple forward trajectory systems (such as those 
arising in implicit integrators, or in model predictive control), as well as, more generally, 
in relations \17\ . 

6 Conclusions 

The breakup of invariant circles in invertible maps of the plane possess a relatively well- 
studied dynamic phenomenology; homoclinic tangles of the invariant manifolds of nearby 
periodic saddle points, which arise from lockings on or off the invariant circle, play an 
important role in the two-parameter picture. Invariant circles for noninvertible maps of the 
plane possess, in some sense, an even "richer" phenomenology. Initially, close to the Hopf 
bifurcation giving rise to them, they are locally invertible, and the Arnold horn scenario 
holds. Large amplitude noninvertible invariant circles, however, may interact with the 
critical curves of the map. This can (and does) lead to new transitions and bifurcations. 
We demonstrated that the one-parameter transition discovered in [16j actually constitutes 
a part of the interior structure of the Arnold horns in the noninvertible two-parameter 
setting; we placed it in the appropriate context by pointing out the role of nearby periodic 
points and their stable and unstable manifolds. The picture we have presented here is 
far from complete, however we hope that our work is a step toward a more complete 
understanding of the transition from smooth invariant circle to chaotic attractor in the 
presence of noninvertibility. 
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